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Abstract: 

A nonperturbative method is proposed for the approximative solution of the exact evolution 
equation which describes the scale dependence of the effective average action. It consists of a 
combination of exact evolution equations for independent couplings with renormalization group 
improved one loop expressions of secondary couplings. Our method is illustrated by an example: 
We compute the /3- function of the quartic coupling A of an 0{N) symmetric scalar field theory 
to order as well as the anomalous dimension to order A^ using only one loop expressions and 
find agreement with the two loop perturbation theory. We also treat the case of very strong 
coupling and confirm the existence of a " triviality bound" . 



1 Introduction 



Nonperturbative field theoretical methods formulated in continuous space and allowing practical 
calculations are relatively rare. Among the best known figure the Schwinger-Dyson equations 
||^] wich constitute a system of exact equations for the IPX vertices. The momentum integrals 
appearing in these equations contain, however, not only IPX vertices but also "bare" couplings. 
Their solution requires an understanding of physics at very different length scales, englobing 
the "short-distance physics" and the "long-distance physics" whose scale is characterized by the 
external momenta of the IPX green functions. The necessity to integrate over a large momentum 
range (within which the relevant physics may vary) makes it often hard to find approximate 
methods which go beyond perturbation theory in some small dimensionless coupling. 
The block spin approach Q in lattice theories seems to be an ideal solution for this type of 
problem: Only the understanding of the physics within a small momentum range is needed in 
order to make the transition from one length scale to a larger one. The effective average action 
Pfc formulates the block spin concept in continuous space. The average action is expressed 

as a functional integral with a constraint which ensures that only quantum fluctuations with 
momenta > fc^ are effectively included. Various field theoretical methods can be used for 
an approximate solution of this functional integral, as in particular, a renormalization group 
improved saddle point approximation This method is nonperturbative in the sense that 

no small dimensionless coupling is required. Nonperturbative phenomena in two and three 
dimensionsal scalar theories [Q], as well as in four dimensional scalar theories at nonvanishing 
temperature @ have been described successfully by this method. 

An infinitesimal analogue of the Schwinger-Dyson equations for the effective average action 
describes the change of Vk with varying scale k. This ensures that only a small momentum 
range ~ /c^ must be contolled. The change of the fc-dependent IPX vertices characterizing P^ 
can now be expressed in terms of such vertices only. Indeed, an exact evolution equation for 
the scale dependence of the effective average action has been proposed recently [0. Xt reads for 
scalar fields tp 

|r,[^] = iTr|(rf [v.] . (1.1) 

Xlere t = In where is a suitable short distance scale and P^^^ is the exact k dependent 
inverse propagator as given by the second functional derivative of P^ with respect to the fields 
if. The trace indicates a summation over internal indices and a momentum integration [|. The 
function contains details of the averaging procedure and causes the integral to be infrared 
convergent. Ultraviolet finiteness is guaranted by an appropriate behaviour of We choose 



1-/1(9) 



- 7271^ (l-^i 



with 

fk{l) = ^^V-^- (1-3) 

We note that for <C A;^ the infrared cutoff ~ Zj.k'^ acts like an additional mass term whereas 
its influence on modes with q^ ^ k'^ is suppressed exponentially. Although eq. (1.1) has the form 



of a renormalization group improved one loop equation (and has actually flrst been proposed in 
this context Q) it involves no approximation. Xt has a simple graphical representation (flg.l) as 

We use 

Tr = Y. J2a = (2vr)"''J7 J d'^qj^a ^^^^^ ^ volume of space-time. For convenience we often 

use Fj. for the second functional derivative devided by fi. 



1 



a loop expression in terms of IPI-vertices only. Equation (^]^) can be related to earlier versions 
of "exact renormalization group equations" by an appropriate Legendre transform |Q. For 
A; — > all quantum fluctuations are "integrated out" and the effective average action becomes the 
generating functional for the IPX green functions (the usual "effective action") in this limit. This 
distinguishes from the "effective action with variable ultraviolet cutoff' whose dependence 
on the cutoff is described by the "exact renormalization group equation" B and which becomes 



a relatively complicated object in the limit /c — > [10|. The simple form of T^, which can be 
described, for example, in terms of a scalar potential, kinetic term and higher derivative terms, 
makes it a reasonable starting point for the development of a nonperturbative method. 
A nonperturbative exact evolution equation (^]^) is not yet a nonperturbative method. In order 
to extract physics one needs in addition a systematic way of solving this equation for /c — > 0, 
with initial conditions specified by the "bare action" Ta at some short distance scale A~^. As it 
stands, eq. (1.1) is a partial differential equation for a function T^ depending on infinitely many 
variables - for example the Fourier modes (p{q) of the scalar field and k. This is, in general, 
impossible to solve. Equivalently, we may expand F^ in terms of invariants with respect to 



the symmetries of the theory. Eq. (O) is then transformed into an infinite system of coupled 



nonlinear differential equations for infintely many couplings (for example IPX Green functions). 
Without a systematic approximation method which reduces this infinite system to a finite system 
(which can then be solved by numerical or analytical methods) not much practical progress is 
made compared to the simple renormalization group improvement discussed in ref. Q. From 
a more systematic point of view a nonperturbative method should be able to compute physical 
quantities and to give an estimate of the error even for situations where perturbation theory 
fails. Our paper contains a proposal for such a systematic nonperturbative method. 
The basic problem we are confronted with can be seen by taking functional derivatives of eq. 



(2) 

(1.1) with respect to 99. We need F^ on the r.h.s. and may obtain an evolution equation for 



this quantity by taking the second functional derivative of eq. (IT). The corresponding fiow 

equation expresses — g|— in terms of F^'^^ , F^^'* and F^."^^ . This feature proliferates to higher Green 
functions: The evolution equation for the n-point function F^"^ involves F^"'"''^^ and F^"'"'''^^ and 

(n) 

the system never closes. For n sufficiently large the couplings F^ are "irrelevant couplings". 
They are usually small in a perturbative context since they are proportional to appropriate 
powers of a small coupling. For nonperturbative problems no such argument is available and we 
have to find a way to estimate F^""*"^^ and r^"^^-* in order to solve the evolution equation for F^."'^ 
The main point for a proposal for a systematic nonperturbative method for scalar field theories 

(n) 

consists in the calculation of F^ for n sufficiently large by a renormalization group improved 

gp(") 

one loop approximation. Xn contrast to the exact expression for which can be derived from 



(1.1) the one loop expression for F^"^ will be approximative. Xt only involves IPX vertices F^™^ 



with m < n. For high enough n the momentum integral in the loop is both ultraviolet and 
infrared finite and dominated by a small momentum range q"^ ~ k"^ R. We therefore expect this 
approximation to be quite accurate if the involved vertices F^"^ do not depend to strongly on 
k. Details of our method and a way of estimating errors will be indicated in the next section. 
We also want to test our proposal for a new nonperturbative method for some problem where 
exact analytical results are known. For this reason we compute in the present paper the /3- 
function for the quartic coupling A of a 0(A^)-symmetric scalar theory in four dimensions. Xn 
section 3 we derive in a short way the general equations. Xn sections 4 to 6 we calculate the beta 
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For smaller n this holds only for J , not for the one loop expression of F^" itself. 
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function of the quartic coupling step by step increasing the number of included couplings. We 
keep all terms in order and find agreement with the standard two loop result. We emphasize, 
however, that we never perform a two loop calculation but only solve the exact evolution equation 
( |1.1| ). We also do not employ here an iterative solution which would be equivalent to a loop 
expansion: All momentum integrals appearing in our calculation have the form of one loop 
integrals. In our approach the two loop result obtaines from an improved one loop calculation! 
Although this result is remarkable by itself it should be seen here as an illustration and check of 
our nonperturbative method. The real power of this method will appear once applied to truly 
nonperturbative problems. A precision calculation of critical exponents in the three dimensional 
theory is in progress and results will be reported elsewhere. As discussed in more detail in the 
next section the proposed method has its limitations for a very fast running of the couplings. 
We will deal explicitly with such a situation in sect. 7, where we discuss the evolution equations 
for very large values of the quartic scalar coupling. The results of sect. 7 confirm the existence 
of a triviality bound in the context of continuum field theory. 



2 A nonperturbative method for scalar field theories 



The evolution equation ( |1.1| ) has a close resemblence with a one loop expression: 



'-Tr 
2 
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—Tr 
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Let us define the " renormalization group improved one loop contribution" T^^^ by 



irrlnl£±^ 

^ In Det (r^') +Rk)-\ In Det {vf + R^ 



and identify Fa with the "classical action" S. Then we can write 

Tk = S + Vf^ + ATfc 



(2.2) 



(2.3) 



where the "correction term" AF^ vanishes in the limit where the last term in eq. (|2.1| ) propor- 

r) (2) (L) 

tional to can be neglected. We emphasize that is already an ultraviolet regulated 

one loop expression and vanishes for k ^ A. This new form of an ultraviolet regulator simply 
subtracts the contribution from quantum fluctuations with momenta > A^ in a way very 
similar to the infrared regularization discussed above. It can be easily generalized beyond the 
scalar theory discussed in the present paper, for example to chiral fermions [11| or gauge theories 
1 12 1, and may therefore be useful for standard one loop calculations in a wider context than the 
present paper. The main advantage of this regulator is the conservation of all symmetries which 
are conserved by the quadratic form 
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ASk = -nJ2ip+{q)Rk{qMq) 



(2.4) 



which provides the infrared cutoff in a functional integral representation of F,fc . 

The disadvantage of the immediate use of fI^^ is linked to the fact that for most theories fI^'' 
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depends strongly on A for ^ ^ oo. This means that the "ultraviolet properties" encoded in 

have a strong influence on r^^\ in contradiction to our strategy to describe physics at 
the scale k only by properties of Tk- For this reason the renormalization group improved one 
loop approximation was used in earlier work Q only for rather than for F^, since ^F^ 

becomes independent of A for A ^ oo. The situation improves greatly, however, if instead of 
Ffc we consider derivatives of F^ with respect to the fields, i.e. appropriate IPI Green functions 
with n external fields. For n sufficiently large the corresponding one loop expression r^"^^^^ 
becomes independent of A in the limit ^ — > cjo ^ and the momentum integral (|2.2|) is completely 
dominated by a narrow range k, k"^ . For k small compared to A we can then take the limit 
A — > oo in order to obtain an expression involving only F^. From the evolution equation 

dt ^ dt 5ipi ■ ■ ■ 5ipn 



r<"' = rf""'' + Ar<"' 
rM(.) ^ irJn-'— i:_i„(rf + (2.6) 



we obtain 



2 1^ 6(pi • • • 6ipr, 

(n) (n) 

We observe that for n sufficiently high F]^ corresponds to an irrelevant operator such that Fj^ 

can be neglected compared to F^"^ for ^ ^ 00. The ultraviolet finitness (for A 00) of the 

momentum integrals in eq. ( |2.6D serves as a direct test for wich F^"^ this procedure of taking 
A ^ 00 is self-consistent. 

Our strategy for a precision estimate of F^ is composed of two parts. First we compute a set 
of exact evolution equations for a set of IPX Green functions F^ by taking the appropriate 
functional derivatives of eq. (^]^). This set denoted by {A} must contain at least all relevant or 
marginal couplings of the theory but can, in general, be larger. We will call these couplings the 
" independent" couplings since they are all treated on equal footing for a solution of the evolution 
equations. The r.h.s. of these evolution equations will depend not only on F^ but also on some 
additional irrelevant couplings F], . As a second step the "secondary" couplings in the set {B} 
are computed by the renormalization group improved one loop approximation (|2.6|). Here it is 



important that F^.'^^ on the r.h.s. of eq. ( |2.6D is computed for a truncation of F^ corresponding 



k 

to a subset of the IPX vertices contained in {^4}. This subset always contains the relevant and 
marginal couplings. More generally, possible choices of subsets of {A} are determined by the 
requirement of ultraviolet finiteness of eq. ( |2.6D . This approximation of F], is then inserted in 
the evolution equation for rj^'^K The system is thereby closed and the evolution equation ( |1.1| ) 
reduced to a finite set of coupled nonlinear differential equations for the independent couplings 
F^^^ which can be solved by numerical or analytical methods. 

Strictly speaking this holds sometimes only for an approximation of ' as discussed below. 
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Furthermore, we may consider r^^-'^^^^ as the zeroth order of an iterative procedure for the 



r D\ 

determination of the secondary couphngs .In fact, we may compute 



approximately by inserting on the r.h.s. the couphngs and Comparing this ex- 



pression with ^r^"'''^^^ computed from eq. (2.6) gives an estimate of the error for the /^-function 



for couphngs . Such an inherent estimate of the error is particularly important for nonper- 
tubative applications of the exact evolution equations where it is not easy to check the accuracy 
of the results otherwise. On the other hand the improved /^-function ^F^"'"''^'' + ^AF^ can be 
used for an even more precise estimate of F^ by solving the evolution equations for the combinet 
set f]."^-'^, F^^^. Further checks and improvements can be obtained by systematically enlarging 

the set of independent couplings F^ for which exact evolution equations are used. 
An obvious limitation of our method arises for situations where the contribution from the term 
■§1^^^^ in eq. ( ^.5| ) becomes comparable or even larger than the corresponding contribution from 
-§lRk- The use of renormalization group improved one loop results for the secondary couphngs 
gives in this case not the correct evolution equations for these couplings. To get an idea when 
this happens in scalar field theories we write F^.'^^ = Z^q"^ + cX{k)ip'^ + . . .. The first term gives 
a contribution proportional to the anomalous dimension. The second contributes substantially 
only for a fast running of X{k). This typically happens for very strong A. We will discuss the 
/3-function for A for the case of strong quartic coupling in sect. 7 and sketch there alternative 
estimates for the secondary couplings. 



3 Exact evolution equation for the average potential 

We consider an 0(A^)-symmetric scalar field theory in arbitrary dimension d. We expand the 
effective average action in terms of the average potential, kinetic term and higher derivatives. 

+ ^5>:n(p,-a2):a^/,+ ...| (3.1) 

where p = ^9?"¥?a and normal ordering indicates that the derivative only acts to the right. Eq. 
(3.1) is the most general effective action which is 0(A^)-symmetric and contains no more than 
two derivatives. The exact evolution equation for the effective average potential can be obtained 
from ( |1.1| ) by expanding around a constant background field 0. One finds 

d , , I f d'^q d , , f N - I 1 \ , , 

W,VM = -Jj:^,j,RM{-ji^ + jf) (3.2) 



with 



Mo(p,(?2) = Zk{p,q^)q^ +Rk{q)+Ul{p) 

Mi{p,q^) = Zk{p,q^)q' + Rk{q)+UUp)+2pUUp) (3.3) 



The neglected higher derivatives do not contribute to 



at 
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and 

= Zk{p,q^)+pYk{p,q^). (3.4) 

Equation ( |3.2[) is a nonlinear partial differential equation for a function of two variables t and p. 
One needs, however, the wave function renormalizations Zk{p,q^) and Zk{p,q^) which must be 
estimated, in turn, by a solution of evolution equations or by alternative methods. We are mainly 
interested in the behaviour near the phase transition, i.e. the massless theory or the theory with 
scalar mass much smaller than the cutoff scale A. Then the evolution of C/fc is mainly described 
@]i by the regime with spontaneous symmetry breaking where the minimum of the potential 
occurs for po > with a k dependent location of the minimum po{k) determined by Ul^po) = 0. 
(Primes denote derivatives with respect to p.) It is convenient to expand the potential in powers 
of p — pq. This transforms eq. ( |3.2| ) into an infinite system of ordinary differential equations for 
the infinitely many couplings {po). In particular, the first and second derivative of eq. ( |3.2| ) 
with respect to p yields 

^Jlip) = T;rf(7V-l)fc'^-2f^If(p)Lf,o(^i)+^;d(iV-l)A:'^4(p)ito'K) 

(3t/^(p) + 2C/f (p)p) Lt^,{w2) + v,k''Z',{p)L^+,\w2), (3.5) 

^C/^(p) = -v,{N - l)k''-^ {K{p)? Li,{wi) - Vdk''-' {W'l{p) + 2U^'{p)pfLi2{w2) 

-2v,{N - l)k''-^Ui:ip)Z',{p)Li%\w,) - 2v,k'-' {WM + 2U^'{p)p) Z',{p)Ltf {w2) 

-v,{N - l)k''{{Z',{p)f)L%\w,) - vak''{{z',{p)f)Lt%\w2) 

+v,{N - l)k'^-'U'^'{p)Ll,{w,) + v,k'^~' (5C/r(p) + 2C/f )(/>)/>) lI,{w2) 

+Vd{N - l)k'Z'ap)Lif{w,) + vak''Z'ap)Lt?{w2). (3.6) 



Here 

and we have used the variables 



= 2'^+^TT2T{d/2) (3.7) 



Wl = Ul{p), 

W2 = U'^{p) + 2U';{p)p. (3.8) 



In terms of the shorthands 



P = Zk{p,x)x + Rk{x), 

P = Zk{p,x)x + Rk{x) (3.9) 
the dimensionless momentum integrals Lf^_^ given by 

^ni,n2(^l,«^2) = A:2("i+"2)-rfy^ rf^^X 2 ±1 | (P + u;i )""! (P + t^s)""' } (3.10) 

where d/dt acts on the r.h.s. only on Rf^ and x = q^. The functions Z'j^{p) and ((Z^(p))^) 
correspond to appropriate moments 

/ dxx2-^Z'i^{p, x)| |(P + ii;i)-"i(P + W2y^} 

i^'^(pK,n. = ° CO , , Z (3.11) 

/ dxx2-^§|(P + wi)-"i(P + u;2)""2| 
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OO , n f 1 

/ dxx2-^ {Z'^{p, x)f I {(P + {P + u)2)-"^} 

{{Z'M)Ynun, = 5S : ^ (3.12) 

/ dxx2~^§i \ {P + wi)-"! (P + wa)-"^ \ 

and similar for Z'i^{p), Zk{p) etc. We have omitted for simplicity of notation the labels ni,n2 
and d as well as the dependence of these moments on wi and W2- Those can be infered easily 
from accompanying functions L^^ ^^{wi,W2). 
The minimum po{k) of the potential is defined by the condition 

U'M = 0. (3.13) 

We are interested in the running of the minimum pQ. For notational simplicity we denote 
u'^\pq), Z^^\pq), Z^f^\pQ) by respectively and understood that Z'^ means the 

moment (|3.12D evaluated at p = p^. Taking the total derivative of (|3.13| ) we get with (^) 

- d du'M -_, 
^=dtP' = — 

= -Vd{N - l)k''-H%{Q) - Vd{N - l)k''Lif{Q)Z'/\ 

-Vdk''-^ (3 + 2C/'Vo/A) L^,i(2Apo) - Vdk''L^f{2\po)Z'/\. (3.14) 
We also define the quartic coupling A = U'^{po) and obtain 

|a = |f/if(po)+^r(Po)^- 

= -Vd{N - l)k''~^\^L%{Q) - vdk''-^ (3A + 2C/"Vo)' 42(2A/Oo) 

-2vd{N - l)fc'^-2^Z'Lg2(o) - 2vdk'^-'^ (3A + 2U"' po) Z' L'^+^ {2~Xpq) 
-Vd{N - l)k''Z'^L%\Q) - Vdk''Z'^Li%\2\p^) 
+Vd{N - l)k''-^U"'L%{<d) + Vdk''-^ [W" + 2t/Wpo) 4i(2A/.o) 
+Vd{N - l)k^Z"Li+\0) + Vdk^Z"L^+,'{2Xpo) 

+U"'6. (3.15) 

This procedure can be extended straight forward to higher couplings. Now it is useful to set 
Zk = Zfc(PO) 0) = Zk in eq. ( |1.2| ) in order to obtain an effective inverse propagator (for massless 
fields) 

Zkipoik), 0)q^ + Rk{q) = J^^- (3-16) 

To separate the running due to dimensions and wave function renormalization from the running 
due to interaction we switch to the dimensionless renormalized couplings (Z = Z^) 



U3 



K = k'^-'^Zpo . 


, A = 


= k'^ 




, M4 


= k 


= k'^-^z-^z' . 


, Z2 


= t 


= k'^-^Z'^Z' . 


, Z2 


= t 



-^z-^x, 

3d-8 2— 4f/(4) 



Zl = k^-^Z-^Z' , Z2 = k^'^-^Z'^Z", 

zi = kr-Z—Z' , Z2 = k^'^-^Z-'^Z" (3.17) 
and obtain the /3-functions 

P.^^ = {2-d-7])K + 2vd{N-l){lf + lf+hi/X) 

+2vd{3 + 2u3K/X)lfsfi2XK) + 2vdli'^^sf+'^{2XK)zi/X (3.18) 
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and 



dX 

/3a = = {d-4 + 27])X + 2vd{N -l)X'^l^ + 2vd{3X + 2u3n)^4si{2XK) 

-2vd{N - l)u3lf - 2vd (5u3 + 2u4k) /f sf (2Ak) 
+Avd{N - l)Xzi4+^ + 4vd (3A + 2u3k) zi4+^s^+^{2Xk) 
-2vd{N - l)z2/f+' - 2vd~Z2lt^^s'i+^{2XK) 
+2vd{N - l)zf4+^ + 2vdzf4+^si+\2XK) 

+U3(^id-2 + rj)K + -^y (3.19) 



Here we have used the expressions 



Ll^{2Xpo)/L%{0) = 4(^) = 4(2Ak). (3.20) 



The anomalous dimension is defined as 



r? = -^lnZfc. (3.21) 



Equation ( 3.19| ) is our starting point for a computation of (3x. So far, no approximation has 



been made and eq. ( |3.19| ) is exact for arbitrary values of A and arbitrary dimension d. In 
addition to A it involves, however, several unknown quantities ^, namely, k, M3, U4, r/, zi, Z2, 
zi, Z2, Zi and z1, and the threshold functions s'^{2Xk). In order to make contact with the usual 
/3-function we have to compute all these quantities as functions of A. We will concentrate on 
the massless theory which corresponds to the (approximate) scaling solution of the evolution 
equation. This scaling solution is characterized by fixpoints for all couplings except for the 
running of A for d = 4. (This running is very slow for sufficiently small A. The true fixpoint is 
the Gaussian for A = 0.) For one of the couplings, which we may associate with /t, this fixpoint 
is infrared unstable. This corresponds to the relevant parameter of the theory, i.e. the scalar 
mass term. For all other couplings (the irrelevant couplings) the fixpoint is infrared stable and 
will be approached rapidly. Due to the existence of the marginal parameter A for d = 4 the 
fixpoints for the various couplings depend on A. Inserting the A-dependent fixpoints for k, 77 



etc. into eq. ( 3.19 ) will then yield fix as a function of A. Using the nonperturbative proposal of 
the last section we compute k, U3 and r/ from corresponding exact evolution equations - these 
couplings plus A correspond to the set {A}. The remaining parameters zi, Z2, zi, Z2, zi and zf 
constitute the set {i?} and are computed using the renormalization group improved one loop 
approximation. The latter can be evaluated as functions of A and k in the approximation of a 
quartic potential with uniform Z}^ = Z^. 

In the present paper we want to check our method for four dimensions where the scalar coupling 
A is small for A; <C A due to the trivality of the (/?^-theory. We want to compute the beta function 



of the quartic coupling A to order O(A^), i.e. we want to know the r.h.s. of eq. ( ^.19 ) to order 



0{X^). Therefore we have to know the expansion of the following quantities in powers of A to 
the specified order 



0(AO) , n3~0(A3), 



^Remember that 2^ stands for an appropriate moment and is in general not equal to the square of z\. 
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Z1~0(A2) , Z2~0(A3), 

~zi ~ 0(A2) , ~Z2 ~ 0(A3), 

r] ~ 0(A2) . (3.22) 

Prom the connection between exact evolution equations and one loop expressions we recognize 
that the orders of the various couplings wich are specified in eq. ( |3.22| ) are also the lowest 
orders which appear in their perturbative expansion. To see this, remember that e.g. in leading 
order of A and at one loop U'j!'{po) is nothing but the contribution to the six-point-function 
at zero external momenta in the symmetric phase and that Z^(po) is a contribution to the q^- 
dependence of the 1-loop four-point-function. Similar arguments hold for the other quantities. 
We see that higher derivatives of the potential (like the (/3®-coupling ^4) or of the wave function 
renormalization don't contribute to the beta function of the quartic coupling to order O(A^). 
Consequently these quantities can be truncated. From the combination of evolution equations 
with one loop expressions we get insight in the important couplings and an efficient tool for 
practical calculations. With this knowledge we have to compute the r.h.s. of 

^ ,(N-l)uslf--^U3ltsf{2XK) 



+^{N - l)Xzill + -^XziIIsI{2Xk) 
^ ,{N-l)z2lt--^~Z2llsl{2XK) 



167r2^ 167r2- 
+2u3K + 0{X'^) (3.23) 

where the fixed point values have to be inserted for all quantities except A. So we have to 
calculate the /3-function of the coupling M3 and the one loop expressions for the wave function 
renormalization constants Zi, Z2, Z\ and Z^- Our aim is to reproduce the known result of the 
two loop perturbation theory, namely 

/3a = Y^(iV + 8)A2 - (^(9iV + 42)Al (3.24) 

For small values of A several simplifications occur. The most important one is that all couplings 
contributing to (i\ ( |3.23| ) in order ~ A^ need only be computed in lowest order in A and that 
hence these corrections are simply additive. This means that for the computation of a given 
coupling, say U3, we can put z\ etc. to zero. This additivity is not crucial for a numerical 
solution of the nonperturbative evolution equations but it helps considerabley for the analytical 
study of the present paper. It allows us to treat the different "corrections" to /3a seperately 
and to compare their relative importance. This will be done in the next sections. A second 
simplification uses the fact that the difference between P and P is a correction of order ~ A^ 
( 3.22| ). In the order we are interested in we can always use ( |3.16| ) 



P[x) = P{x) = (3.25) 

in the integrals ( |3.1C1| ) and omit the /c-dependence of Z^, in ^ . This permits a simple analytical 
computation of the constants [^. In particular, one finds 



l\ = l, ll = l, ll = 3ln(J\, If = 2, /6 = | (3.26) 
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Finally a third simplification concerns the threshold functions s^(2Ak). They can be expanded 
for small Xk according to 

lisi{2\K) = li-2nXKli^, + ... (3.27) 
We observe that we can put sf^ = 1 except for the term ~ A^s^ in eq. ( |3.23| ). 



4 The average potential with uniform wave function renormal- 
ization 

We first consider in this section a uniform wave function renormalization ^fc(P) Q^) = ^i^^p, q^) = 
Zk and therefore set 



Z' = Z" : 

Y = Y' 



0, 

: 0. 



In this approximation one obtains from differentiating eq. ( p.6D with respect to p 
d 



(4.1) 



(4.2) 



The lowest order contribution for the running of k and u-^ therefore read 

(2 - d)K + 2vd{N + 2)/f , 
{2d - 6)^3 - Avd{N + 26)/|A^ 



dt 



d_ 
dt 



Us 



(4.3) 



The fixed points are now easily calculated as the roots of eq. (O). One finds in leading order 

N + 2 



= 2Vd 



d-2 



■n 



which is infrared unstable and 



2vd 



Note that the last value coincidences with the one loop expression for 



u 



(1) 



2{N + 2&)Vdy / dxx^-^p-\x) 



(4.4) 
(4.5) 

(4.6) 



which follows from eq. (|2.6D at zero field = and n = 6. This coincidence can be understood 
by looking at the estimate of error (2.7) of the /3-function caused by the use of "renormalization 



group improved" one loop expression ( [4.6D instead of the evolution equation (O). Comparing 
the estimate of error for us 



^ = ^I^T^h^ = 01 = \r 
dt dt ^ ^ 2 (5(^i(g = 0)---<5(/?6(g = 0) 



with its exact evolution equation 



dus 



d 



^ = ^ri% = o] = -T. 



5(pi{q = 0) • • • 5ipQ{q = 0) 



1/3 = 

(4.7) 



(4.8) 



1^=0 
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we recognize that the estimated error is of higher order in A since ^r^^^[0] 



whereas 



We also observe that the difference between rj^ [99 = 0] and T)^' evaluated at the 



.(6) 



minimum is of higher order in A. It is therefore no surprise that the lowest order contributions 
eqs. (4.5), 
dimensions 



|) coincide. Inserting the fixed points (4.4), 
1 „.,o 1 



in eq. (3.23) we obtain in four 



Px 



167r2 



(iV + 8)A' 



lf4 (20A^ 



{Wit 



2\2 



which is numerically 



/9a 



1 



167r2 

to be compared with eq. (3.24). 



(A^ + 8)A^ 



(167r' 



2^2 



(17.26iV + 75.95) A^ 



(4.9) 



(4.10) 



5 Field and momentum dependence of the wave function renor- 
malization 

We next take the field and momentum dependence of the wave function renormalization into 
account. We want to compute zi, zi, Z2 and Z2 in a renormalization group improved one loop 
approximation. For this purpose we write Zk{p,Q'^) — Zk{0,0) as an appropriate difference of 
two point functions and similar for Z^.. This difference is evaluated according to ( p. 61 ) where on 
the r.h.s. an approximation for T^^^ is infered from the ansatz (3.1) with Zk{p, —d^) = and 
Yk{p, —d'^) = 0. In a first approach we neglect the momentum dependence and approximate 
^kip) by Zk{p, q'^ = 0). The one loop expressions are derived in appendix A and read 

Zfc(p,0) = Zk{0,0)+8^k''-'{U'^{p)fpEl^{wuW2), 

Zk{p,0) = Zk{0,0)+4^k''~'{3Ui:{p) + 2U'^'{p)p)^pE^4wi,W2) 

+4{N-l)^k'^-'{U'^{p))'pEl,{w,,W2). (5.1) 
Here we have introduced the dimensionless integrals 

Ei,,n,{wi,W2) = k^^''^+''-'-^^"^ dxx2P\P + w,)-''^{P + W2)-''' (5.2) 

^ 

where wi^2 are given by ( |3.8| ) and we have used x = q^, P = ^ with P given by ( 3.25| ). The 
p-derivatives of eqs. (5T) at fixed po yield in lowest order 

Z' = S^k'^-^X^Ef, 



Z" = -US^k^'-^X^El 
Z' = 4{N + 8)^k'^-^X'^E: 



- -4, 

Z" = -32{N + 2Q)^k'^-^X^El (5.3) 

with E'^_^_^_^^ = E!^_^ n2i^^ 0). We switch to the dimensionless renormalized quantities and find for 
d = 4 

_ 1 \2 4 

- 16^^'^' 
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Z2 



N + 8 



327r2 
iV + 26 
4^ 



where 



Z2—n d 
En 



and 



4 



In 



(5.4) 
(5.5) 



1 

2' 5 2 

Again these one loop expressions are in lowest order equal to the appropriate fixpoints of the 
corresponding evolution equations. As before this holds because we include only the leading 
order contributions to the one loop expressions or the appropriate evolution equations. 
The formulae ( |5.4D are only approximations and do not include the full contributions in leading 
order in A. For a precision estimate in this order we have to use approximately weighted functions 
as indicated by eqs. ( p.ll|) , (|3.12| ). We therefore have to take the momentum dependence of 
Zf^{p,q^) into account and employ the renormalization group improved one loop expressions 
([A. 61) given in appendix A. We notice that ^''^'^^f^ ^ = 1 + O(A^). This is the well known result 
from standard perturbation theory. In the /?A-function (|3.23| ) the quantities Z', Z' enter up to 
order O(A^) and Z",Z" enter up to order O(A^). Deriving equations ( [A.6| ) with respect to p at 
the minimum po we obtain in leading order 



4'(po,g') 



-2|(2.)- 



16|(2vr)- 
-(iV + 8)^(2vr)-'^ 



d'^qP-\q) [p-\q + Q) - P-\q) 
d'^qP-\q) [p-\q + Q) - P'Hq) 

(fqP-^{q) \p-^{q + Q) - P'^iq) 



A{N + 26)^(2^)-^ / dUP~\q) \p-\q + Q) - P~\q) 



(5.6) 



We notice that in the required order in A the mass term ~ 2AkA;^ in the denominators can be 
neglected and we have only to solve the one loop integrals for a massless theory. We evaluate 
these integrals using the Laplace transform 



and obtain 



p-\q) 



1 



dae 



aq 



(5.7) 



d^qP-Hq) \p-\q + Q)-p-\q) 



dadf3 







1 (a + /3)-2 



(5.8) 



and 



J(Q2) = 1(2^ 



d^qP-\q) [p-\q + Q) - P-\q) 



dadpd'-f 







«+/3+7 -lUa + P + j) 2. 



(5.9) 
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After some transformations we obtain in four dimensions (See app. B) 



and 

j(Q2^=^ k-^y- — til! — I [' dz-^^^ - r dz-^—-^] . (5.11) 
^ ' ^,in + 2)lin + l)\k^ J \Jo (l-z)2 ^ {I - z)^ j ^ ' 

The dimensionless renormalized quantities read in lowest order in A 

oo 



(n + 2)!(n + l) V2 
(_l)n 




~z,iQ^) = (iV + 8MA^E^^^^,^^^^^^^ 
^0 (n + 2)!(n + l) \ 



(-1) 

^ ' ^ ' ^^{n + 2y.{n + l) 

Notice that the (n = 0)-terms of the sum correspond to eqs. (|5.4| ). We use these expressions to 
calculate the appropriate moments ( |3.11| ) which occur in the /9A-function (3.23). 



' ^^{n + 2)\{n + l)\2l 



6+2n 
2 



= (241n3 + 161n2-201n5-4)^;4A^ 

= (iV + 8)(121n3 + 81n2- 101n5-2)u4A^ 

,6+2n f n r. H r. 

)!(n + l) 1 
16(12 In 3 - 5 In 5 - 8 In 2)^;4A^ 



Z2^? = 32t;4A3 y '-^ -/6+2n / dz-^-— - / ' dz- 

' ^Q(n + 2)!(n + l) 1 [io (1 - J I (l - zY 



z^lf = 8{N + 26)t;4A3 V / rfz-^-— - / ' ^ 



„t^^(n + 2)!(n + l) ^ ^0 (1 - (1 - z)^ j 

4(A^ + 26)(12 In 3 - 5 In 5 - 8 In 2)7;4A^. (5.12) 



(See app. C for the evaluation of the sums.) Comparing with the expressions ( ^.41 ) we see that 
the Q^-dependence of the wave function accounts for 21%, 21%, 30% and 30% of the moments 
zi, zi, Z2 and Z2, respectively. Inserting the results ( 5.12 ) into the /3A-function ( p.23| ) reads 

In order to compute the complete coefficient ~ A^ we only need now the anomalous dimension 
r] in order A^. 
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6 The anomalous dimension 



The anomalous dimension rj is defined by identifying with Zi^(po(k), = 0) 



-|lnZ,(po) 



(6.1) 



where 



Zk{p) = Zk{p,Q^ = 0) 



Q-^ Hm 



d 



hm 



d „(2), 



(6.2) 



is evaluated for a constant field configuration. With the exact evolution equation eq. ( |1.1D we 
obtain an exact equation for the /c-dependence of Zk{p) (for iV > 2 Fl) 



dZkip) 



1 



lim 



d 



-Tr 



+ Rk] dtRk 



We perform the functional derivatives (omitting the internal field indices) 



(6.3) 



dt 



1 



lim 



d 



Q2^o 9Q2 6ip^Q)dip*2{Q) 



Tr 



— i7 ^ lim -r—^Tr 

2 Q2^o 



(g,-g) + i?,,((Z) 



r(')(g,-(7)+i?fc(<7)] dtRk{q) 
1 



rr(Q,-Q> ('7,-g) + ^fc('7) '5ti?fc(g) 



+2 



r(')((?,-g)+iifc(g) 'r^^^(Q,g,-g-Q) r^^^(g + Q, _ Q) + i?,((; + Q) 



.(3) 



(2), 



with 



ri')(g,-<7) + i?,.(g)l ^dtRkiq) 



(6.4) 



(6.5) 



(5v9(gi)---599(g„) 

containing a factor ~ (5(gi + . . . + g„). Inserting /> = /fjQ in eq. ( |6.4| ) and combining eq. ( |6.1| ) 
with eq. (|6.4D gives an exact nonperturbative expression for the anomalous dimension r] ! 



It involves the /c-dependent full propagator {t'^^'' + Rk^ , the IPI three point function for 
arbitrary momenta and the IPI four point function for two pairs of opposite momenta. One 
of the momenta (corresponding to Q'^) can be taken small and only an expression quadratic in 
this momentum is needed. (We observe that eq. (6.4) can easily be generalized into an exact 
evolution equation for the full momentum dependent "wave function renormalization" Z^^p, Q'^) 
by omitting the operation limQ2^o subtracting from the now Q dependent r.h.s. the value 
for Q = and deviding by Q^, i.e. by combining eq. (1.1) and eq. ( A. 2 ).) 

The ansatz ( p] ) only describes T^^\o,q, —q) and T^^Uo,0,q, —q) or other combinations with 

®For A'' = 1 one should use functional derivatives with respect to (p = (pi . This changes the formulae of this 
section but leads to the same final result for rj. 
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only two nonvanishing momenta since higher derivative terms of the structure d'^ipd^pdip and 
{dip)^ are not included. In a first step we neglect the momentum dependence of r[,^^ and F^.^"* 
except for the contributions quadratic in the momenta. We can then determine the r.h.s. of 
eq. ( |6.4D from the ansatz (3J) with Zk and Yfc depending only on p. In this limit ry has been 
computed in Q,[13| and we quote only the results. The scale dependence of Z^^p) at fixed p 
reads in this approximation 



dZk{p) 



dt 



dt 



dt 



(6.6) 



with 



dZkip) 



(a) 



dt 



dZk{pf^) 



dt 



Vd [{N - l)Z'k{p) + n(p)] k^-^L%{wi) 
+Vd [Z'.ip) + 2Z^'(p)p] k''-^Lt4w2), 



(6.7) 



+Vd{Yk{p)f pk''-^Qil{wi,W2) - 8vdZi{p)U'^{p)pk^-^Ll^iwuW2 



4fd 



(Z^(p))^pA;'^-2L5'f(Tz;i,«;2)-4t;,Z^(p)n(p)pA:'^-'^?f («^i,tx;2) 



+ 



d 
d 



d-2 Td+2i 



Z'dp)U'ap)pk''-'Nl^{w^,W2) 



+^Z',{p)Yk{p)pk'-'N^+\w,,W2). 



(6.8) 



Note that eq. stems from the part of (|6.4| ) which involves the four point function whereas 

eq. (6^) is originated from the part of ( |6.4| ) which involves the three point functions. The 
integrals are defined by, with P — 



dx 



and the integrals 



d_ 
di 



(6.9) 



d_ 
'dt 



2x 



dxx 2 
4x 



d d 



{P + Wl)-'"^{P + W2) 



are related by partial integration to other integrals 

d 



A^ni,„2(«^i,^2) = A:2("i+'^2-i)-rfy^ dxx2^{p2(p + ^^,)-"i(P + y;2)-"2 



} 



(6.10) 



(6.11) 



through 



2ni 



^ '-M^'^ll^^{wr,W2) + ^M„^+2.^,(u;i,^2) 



+- 



2n2 



2a 



^Ni+^^-\w,,W2). 



(6.12) 
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Note, that the £'^^ -integrals are closely related to the M^^ -integrals since §ig{-§i) = 
—2x-^g{-^) for an arbitrary function depending only on the ratio p-. In lowest order of A 
eq. (B^) reads at the minimum po 



^ = VdiNZ' + Y)k^-^Ll,{0) + ^vA'pok''-'Ml,{0). (6.13) 



Expressed in the dimensionless renormalized quantities (and y = Z '^k'^ "^Y) one obtaines for 
the anomalous dimension ( |6.1| ) 

= iQ^^mfKX^ + 2vd{Nzi + y)lf -zi(^ + {d- 2)k) (6.14) 

where 

mi = -iz"~2M„%(0). (6.15) 
We remember y = zi — zi and write at the fixpoint of our theory 

= le^m^K^A^ + 2vdli{ii)* - 6vdlf{zi)^. (6.16) 

Since we have neglected the momentum dependence of Z^ipjQ"^) the appropriate one loop ex- 
pressions are given by ( |5.4| ) and we obtain in the present approximation and for d = 4 

Next we include in eq. ( |6.4D the momentum dependence of T^^^ and T^j^^ beyond the approxima- 
tion quadratic in momenta. We observe that the momentum dependence of the three point and 
four point function is always of higher order in A compared with the momentum independent 

part. The contribution ~ (r[,^^^ is already ~ A^ and we therefore do not need the momentum 
dependence of the three point function for a computation of rj in order A^. On the other hand 
we need the quantity 



Ak{Q,-Q,q,-q,) = T)^\Q,-Q,q,-q) -r'^^> {0,0, q,-q) 

f\Q,-Q,0,0)+4 



rl^^(Q,-Q,o,o) +ri')(o,o,o,o) (6.18) 



in order A^. This summarizes the neglected terms and the insertion of ( |6.18| ) into ( |6.4| ) gives the 
complete correction to rj in order A^. In fact, the relevant Q-dependence of the first term on the 
r.h.s. of eq. (|6.4|) is given by an expansion of rl"^^ quadratic in Q. The most general form of the 



k 

part of r^^\Q, —Q, q, —q) can be written as 

rl'^(Q, -Q, q, -q)=F + G{q^ + Q^) + H^{q^) + Hr{q^)Q'' + H2{q^){qQf + O(Q^) (6.19) 
with 

i?o(0) = 0, H{0) = 0, Hi{0) = 0. (6.20) 

The first two parts involving F and G have already been included in the approximation used 
in the first half of this section whereas HQ{q^) does not contribute to t]. It is easy to verify 
that the remaining term Hi{q^)Q'^ + H2{q^){qQ)'^ is just the part quadratic in an expansion of 
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^k{Qj ~Qi Qi ~q)- As an illustration of this argument we give the expansion of F^"^^ in fourth 
order in the momenta for = 1 

.(4), 



r):>iQ,-Q,Q,-Q) 



3Ui:{p) + 2pUl!'{p)+4p'ul'\p) 
+ iq' + Q')iZ'k{p) + 2pZ'^{p)) 
+ {q^Q^ + 2{qQf)Wk{p) 



(6.21) 



and note that the last term ~ Wk{p) corresponds to the contribution from neglected so far. 
We will consider as one of the "secondary" couplings in the sense of sect. 2 and evaluate it by 
a " renormalization group improved" one loop expression. In lowest order we can neglect the 
mass terms ~ IXnk^ in the propagators and only need a one loop calculation for the massless 
case, assuming a purely quartic potential and a constant wave function renormalization Z^. 
With these approximations the one loop expression for the four point function reads 



{r'j^\Q,-Q,q, 



22aa 



}?diag{{N + A)A + 2B{q, Q),{N + S)A + {N + 8)B{q, Q), 



{N + 4)A + 2B{q, Q),...,{N + 4)A + 2B{q, Q) 



(6.22) 



Af-2 



where 



A = ^i2n)-'' I d''pP-\p), 



B{q,Q) = -{27^)-'' j d''pP-\p)[p-\p-Q-q)+p-\p-Q + q) 



(6.23) 



The expressions A and B are both ultraviolet divergent in four dimensions but these divergencies 
cancel for the difference ( 6.1^ ) 



Af-2 



Afc((3,-Q,g, -g, 



-X^diag 



1 



2,iV + 8,2,...,2 I ^(2^)-^ 



V 



d'^pP^^ip) 



2P-\p) +p-\p-Q-q) + p-\p -Q + q)- 2P-\p + Q) - 2P-\p + , 



(6.24) 



which is finite in any dimension less than six. Again we evaluate this expression using the 
Laplace transform (|5.7|) and obtain 



^k{Q,-Q,q,-q,] 



-VdT ( ^ ) X'^diag 



/ N-2 \ fc-2fc- 

2,iV + 8,2r^ J J 



dadp 



a/3 ^2 a/3 ^2 

l-2e "+/'^ +2-2e -+'3^ +2 



(« + /?)" 



(6.25) 



With the /((5^)-integral ( 5.10| ) we write in four dimensions 

/ 



Af-2 



Afc((5, -Q, q, -q, ) = -V4X diag 

- / , „ . 9 \ n+l 



A;2 




-1)"+1 



^Jn + 2)!(n + l) 



2, 



9 \ n+l 

q^\ 



(6.26) 
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We now insert this expression into eq. (|6.3|), perform the UmQ2^Q operation under the 
momentum integration corresponding to the trace and obtain with (|6.1|) the contribution 



/•oo 

A77 = -2,{N + 2)vl\^ dydtP{y)P-\y)Y^ 



n+l /T,\n+1 



4)n+l /y 



,n + l)! V2 

71=1 ^ ' ^ 

where we set y = /k'^ and used eq. ( 3.25| ) to perform the integration. Adding this contribution 



Ar] to the result from eq. ( 6.17 ) we obtain in order 

N + 2 A2 



(16^2)2- 



(6.28) 



This coincides with the perturbative two loop result |T^. Insertion of ( |6.28 ) into eq. ( |5.13| ) 
yields for the /?A-function 



_N + 8 , 9iV + 42 3 



Again, the two loop coefficient IQ is exactly reproduced. Our "nested one loop calculation" 



gives in next to leading order in A exactly the same results as a two loop calculation! 

7 Large quartic scalar coupling 

In the preceding sections we have concentrated on small values of A. Our nonperturbative 
exact evolution equations are by no means restricted to this case. We demonstrate this here 
by discussing the leading order contribution to Px for very large values of A. We also combine 
approximately the two regimes for large and small quartic couplings in order to obtain the /3- 
function for the whole range of A. This is justified by the existence of only a small transition 
regime which is reached quickly for scales only somewhat below the ultraviolet cutoff. Our 
starting point are the exact formulae ( |3.18| ) and ( p.l9| ) for /J^ and (3x which we combine to 

Px = {d-4 + 27])X + 2vdiN -1)X'^4 

+2vdiN - 1) [2Xzili+^ + U3Zi/Xlf+^ - Z2/f+' + zllt^} 

+2vdlfsf{2XK) I2U3K/A - 2u3 - 2u4k| 
+2i;rf/f+2^f 2(2A^) {ims/A - Z2} 

+2ud| (3A + 2u3Kf lisi{2XK) + 2 (3A + 2u3k) ziZ^+2^^+2(2Ak) 
+5?/fSf^(2AK)l. (7.1) 



Our strategy for an evaluation of (3x will consist of a series of assumptions on the maximal value 
(in powers of A) that different quantities on the r.h.s. of eq. ( [7.1D can take for very large values 
of A. From there we compute the leading contribution to Px- In a second step we establish 
the self consistency of our assumptions by discussing explictly the evolution equations for the 
quantities of concern. 

Eq. ( |3.18|) gives no indication that k should be suppressed by inverse powers of A and we 
assume in the following k = 0(1) in the sense of counting powers of A~^. Then the threshold 
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functions s^(2Ak) give suppression factors ~ A"^""'"^). We will restrict the discussion of this 
section to > 2 where the contributions ~ A"^"'"'"^^ are small compared to the Goldstone boson 
contributions ~ (A'^ — 1). Neglecting terms ~ s'^{2Xk) in a first approximation we find 

f3x = 2vd{N -l)X^4 + {d-A + 27])\ 

+2vdiN - 1) {2Azi/f 2 + n3Zi/A/f+2 _ ^ ^2;d+4| ^ (7 2) 

We remind that the quantities Z,^ (|3.20| ) have to be evaluated here with P = Zk{p, (p')q^ + Rkio^) 
according to ( |3.9| ). The difference between this and the values for small A (where P Zi^q^ + 
Rk{q)) as well as the last term in ( [7.2| ) can be attributed to the wave function renormalization 
Zfc(/9, Q^)/Zfc. We therefore need an estimate of this quantity. 

Let us start at some ultraviolet cutoff A with Za{p, Q'^) depending only weakly on p such that 
\Zj^{p,Q'^)\ < 0(1) and |Z^(/9, Q^)| < 0(A). We first assume that this situation remains valid 
for all scales k for wich A remains large. Assuming also that |n3| is at most O(A^) we find the 
simple leading order expression 

/3a = 2vd{N - l)A2/f + (d - 4 + 2r?)A. (7.3) 

Here the unspecified Q^-dependence of Zk{p, Q"^) only affects the value of I2 but not the general 
structure of eq. ( |7.3| ). Unless the Q^-dependence of Zk{pQ,Q'^) is extremely dramatic will 
always remain of order one. As long as -^^^I^HQ, —Q,q, —q) does not exceed a value ~ 0(1) 

and g^^r^ \Q,Q, —q — Q) 0(1), ^^T^ {Q + q) ~ 0(1) we conclude from the exact relation 
(|6.4| ) that r] remains of the order 0(1). With all these assumptions one finds the leading order 
contribution 

/3a = 2vdiN - l)X^li (7.4) 

Neglecting in addition the Q^-dependence of Z}^{pQ, Q"^) the coefficients I2 are given in appendix 
C with I2 = 1. Before a discussion of the validity of our assumptions several comments are in 
order at this place: 

i) For large N we recover the result of the leading contribution in the large A^ approximation 



Px = 2vdNX^4 + {d-A)X, 
rj = 0. (7.5) 

We remind that the large A'^ approximation obtains formally by neglecting a collective excitation 
which roughly corresponds to the radial mode. With our assumptions the contributions from 
the radial mode are small because of the threshold functions s^(2Ak) and similar suppressions 
from the propagator of the massive radial mode ~ A~^ in other expressions. This motivates the 
qualitative correctness of ( [7. 4) ) also for moderate values of A^ (except A^ = 1). 

ii) We note that there is no dramatic difference in the behaviour of (3\ for small and large 
A. We may therefore extrapolate for intermediate values of A by connecting (3x for small and 
large A at the transition value Xtr where they take equal values. Roughly speaking Xtr marks 
the transition from the perturbative to the nonperturbative regime. We define 

/3 = (/3a - (d - 4)X)/i2vJiX^) (7.6) 



19 



such that, for d = 4, 

^ = P + 8-tI^A forA<A,. .77. 
I + 7 for A > At/ ^ ■ ^ 

(The following discussion will remain valid with minor quantitative modifications in arbitrary 
dimension d. We also allow for a correction factor 7 accounting for a possible modification of 
the leading order result ( [7. 3D which corresponds to 7 = —1.) From ( |7.7| ) we obtain the transition 
point 

9iV + 42 ^ ^ 

(For = 4, 7 = —1 the transition point is A^^ = 18.2.) At A^r the ratio of the two loop to one 
loop contribution for f3 is given by 

^^ = -1^ (79) 

A^ + 8 ^ ■ ^ 

which is of the order one except for very large N. 



iii) The precise form of /? for large A, e.g. the value of 7 as long as 7 is sufficiently large 
compared to — A^, is not very important for the qualitative behaviour of our model. Taking into 
account the leading contribution 



— = 2vdi{N + ^)\^ (7.10) 



one finds for the scale dependence of A 

For d = 4 the value Xn- for the transition to the perturbative behaviour is reached for k > kt^ 

9N + 42 

(The scale ktr is obtained for A^^(A) = 0. For positive values of A^^(A) the transition occurs at 
higher values of k. This is used below to establish an upper bound on the Higgs scalar mass.) 
For A^ = 4 one finds 

fctr-^r 0.056 for7 = -l 

A 1 0.087 for 7 = ' ^ ' 

Obviously the theory reaches the perturbative regime very fast. In four dimensions there is no 
(/9^-theory with large renormalized quartic coupling at scales much smaller than the ultraviolet 
cutoff. For k below the scale kpt where perturbation theory becomes valid we use the solution 
of ( |6^ with Xpt = \{kpt) 

1 , 9A^ + 42 , (1%tt'^{N + 8) 1 ^ 



X{k) 167r2(A^ + 8) \ DAT + 42 \{k) 



1 9Af + 42 , + 1 \ A^ + 8, A;„i 

1 TT-, T m 1 H TT- m 

Xpt 167r2(A^ + 8) \^ 9Ar + 42 Xpt j IQn^ k 

N + ^K 9A^ + 42 9A^ + 42 , A^ + 7 ^ 

= In 1 In (714) 

167r2 k 167r2(A^ + 7) ^ 167r2(Af + 8) 8-7 ^ ' 
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For the last equation we have inserted kpt = hr ( 7.12| ). This constitutes an upper bound for 



possible values A(fc) for arbitrary scales k < ktr- We observe that this upper bound does not 
depend strongly on 7, especially for large ratios K/k. Let us compare the result ( |7.14 ) with the 
result of high order strong coupling expressions on the lattice [15|, namely 



In -?-rM^) =^7?-r ln- + C , (7.15) 



\{k) 167r2(iV + 8) V 167r2 ' 'J levr^ V k 

where C = 1.9 for = 4. We observe agreement in the large N limit and also in leading order 
~ A^^. The constant C computed from ( [7.14| ) obtains as C = —2.8 and deviates from the lattice 
results. This may due to our neglection of nonleading terms in the /3-function for large A and 
in an inaccurate matching between the strong coupling and weak coupling regimes. Also the 
lattice definition of A is not identical with our definition. If we include in the strong coupling 
/3-function ( [7.4[ ) a nonleading term 5 ■ A and determine 5 such that at \tr both 15\ and match 
for the strong and weak coupling regimes, we find (for 7 = —1) 6 = 0.26, C = —1.8. 

iv) The upper bound on A(A;) translates into an upper bound for the scalar mass both in the 
symmetric phase and in the phase with spontaneous symmetry breaking. This bound depends 
on the ultraviolet cutoff A. In the phase with spontaneous symmetry breaking the value of the 
renormalized field corresponding to the potential minimum goes to a constant pQ larger than 
zero for A; ^ 0. (The scaling solution with almost constant k for small enough A ceases to be 
valid for k smaller than some characteristic scale of the order /c^ ~ X{k)pQ.) The renormalized 
mass is then approximately givenj^ by 

m2 = 2A(/9o)po. (7.16) 

With ~ 246GeV we obtain an upper bound M^ax as a function of A by inserting the 

upper bound for \{pq) from eq. ( [7.14 ). It is qualitatively similar to the lattice results [p^] , [|l^] 



and deviates only in an overall factor of the constant C . 

V.) One should note that (3\ can be written as a function of A on a given trajectory, i.e. 
for given initial values of T\. This is achieved by expressing all other couplings on the r.h.s. 
of ( |7.lD in terms of A by inserting the solution A(A:), i.e. n3(A) = U2,{k{X)) etc., provided A(A;) 
is a monotonic function. We emphasize, however, that this function /3a (A) is not a universal 
function. In general the details of /3a (^) depend on the particular initial values of the couplings 
specifying Fa. Only in the vicinity of a fixpoint where A varies only slowly, k is fixed by its 
(A-dependent) critical value corresponding to the critical trajectory, and all other couplings take 
their values at infrared fixpoints as functions of A, the function /3a (A) becomes universal. In 
four dimensions this happens for small A (the vicinity of the Gaussian fixpoint). In contrast, 
for very large values of A the evolution is extremely rapid and one is far away from any fixpoint 
behaviour. Details of /3a (A) for large A are therefore indeed expected to depend on detailes of 
the short distance action Fa. The problem of estimating /3a (A) for the most general form of Fa 
is certainly not solvable. In extreme cases Fa may correspoond to a theory in a universality 
class completely different from the (/9^-theory. For example, the effective action Fa for a theory 
with light fermions and Yukawa couplings to the scalars can be written in a form of a (non- 
local) action involving only scalar fields by integrating out the fermionic degrees of freedom. 
This would certainly lead to a different function /3a (A) even for small values of A. (Additional 

^We define the mass for nonvanishing momentum in order to avoid infrared problems from the Goldstone 
bosons. Since \{k) changes slowly for small enough A we can take A(po) instead of the appropriate four point 
function at nonvanishing momentum with a good accuracy. 
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contributions from Yukawa couplings would arise.) We will be satisfied to establish that the 
qualitative behaviour of ( |7.7| ) holds for a large class of initial values T\ within the range of 
attraction of the Gaussian fixpoint. The discussion of the validity of the assumptions leading 
to ( |7.7| ) will also give information on this range of attraction. For initial values in the vicin- 
ity of the boundary of the range of attraction of the Gaussian fixpoint we expect qualitatively 
different evolution. No "triviality bound" of the type discussed before can be derived in this case. 



Let us now discuss the assumptions leading to ( [7. 7]) in more detail. The most important 
one was that k is not of the order or smaller such that contributions from the radial mode 
~ s5^(2A/t) can be neglected. For k ~ A^^ the leading contribution to the evolution of k would 
be given by ( ^.18 ) 

"-'^l + il). (7.17) 



dt 167r2 V A. 

This would drive k very quickly to zero and we conclude that values k ~ A^^ or smaller cor- 
respond to the symmetric phase. We are interested in the phase with spontaneous symmetry 
breaking near the phase transition. Here k cannot be smaller than of order one. Our next worry 
is the value of the dimensionless (/3^-coupling U3 . If M3 would be of the order A^ the suppression 
due to the threshold functions is insufficient and (3\ would obtain contributions ~ A'^. The 
validity of the approximation (^^) requires that U3 is at most of the order A^ . At first sight this 
seems to contradict the one loop estimate ([4.6[) . We should remember, however, that large values 
of A correspond to a situation with rapidly varying couplings, in contrast to the approximative 
fixpoint behaviour for small and slowly varying A. For an estimate of U3 we integrate its leading 



order evolution equation which may be obtained from eq. (3.6) after derivation with respect to 



p and introduction of the dimensionless renormalized couplings ( 3.17] ) 



^ = -4vdi{N - 1)A3 + 6va4{N - l)An3. (7.18) 



As long as is small compared to A^ we can use the solution ( 7.11 ). The first term in ( [7.18| ) 



leads to an increase in ns until it reaches a value ~ A^. Then the second term in ( [7.18 ) becomes 
important and u^/\^ is driven to a zero of the r.h.s. of 

a^X3_ 2^n3 3^A _ N-l ( , .U3\ 



(for d = 4, 7 = —1). We note that the fixpoint 



corresponds now to a value of of the order A^. The difference compared to (|4.5| ) arises from 
the second term in ( [7. 18) ) which has been neglected in sect. 4. By similar arguments it is easy 
to establish that U4 is at most of the order A^ in agreement with our approximations. 
Let us finally make an estimate for the moments and Z2. The one loop estimate ( 5.12| ) suggests 



values zi ~ A^, ^2 ~ A^ in contradiction to the validity of our approximations. If we write the 
evolution equations for Zi in the form 

^Zi = -A{X)+Bi{X)zi + ... (7.21) 
ot 
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only the terms Ai ~ ■j^,A2 ~ are directly related to the one loop formulae of sect. 5. 
(For the discussion of section 5 it is implicitely assumed that Bi is dominated by the canonical 
dimension of Z', Z" . This is only valid for small A.) In leading order the evolution equation for 
4(po,Q' = 0) and Z^'(po,Q' = 0) Q, (|D yields 



(ziK) = 2vd\ (^mi + {N- l)4zif^ , (7.22) 



di 
d 

— (z2k) = 4Wrf(iV - l)/fz2«A 



-^vA^ {^mi + {N- l)liziK 



+2vdU^ Qmf + {N- l)4zit^ . (7.23) 

Here we have approximated L'[^{wi,W2) — ■^L'[q{wi) and similar for other threshold functions 
and we have omitted terms ~ 1^. (This becomes exact in the large limit). We conclude that 
for large A there are infrared fixpoints in zik and the ratio Z2k/A 

A A {N-l)dli\ 4 iV-lV 3 4/ ^ ^ 



(Here we have inserted the fixpoint value for ( 7.20 ) and the last part of the equations is 
evaluated for d = 4 with constants l'^,rn^ given in appendix C.) From this we infer that \zi\ 
grows within the regime where A is large at most to a value ~ 0(1) and \z2\ becomes at most 
~ 0(A). In consequence the last term in ( |7.2| ) can indeed be neglected in leading order in the 
regime of very strong quartic scalar coupling. With the estimations 

|k|<0(1) , \u^\<0{\^) , \ui\<0{\^), 

\zi\<0{l) , |z2|<0(A) (7.26) 

and s5^(2Ak) ^ 0(A^^"^^^) we find that the approximation ( |7.3D is indeed valid provided 

\zi\ < 0(A2) , \z2\ < 0(A3). (7.27) 

We therefore quote the evolution equation for Zk wich has been computed in sect. 5 of 0] and 
reads in leading order for large A 

-2vdk''~\N - l)Yk{p)U'^{p)pL%{U'^{p)). (7.28) 

Expressed in dimensionless renormalized couplings we obtain for the scale dependence of Z and 
its derivatives in leading order [zq = Z /Z) 

+Avd{N - l)liyoXK, (7.29) 
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= _2v,{N-l)lfy2K-^{N-l)miusXK 
+4vd{N - l)4yu3K + GvdiN - l)/fyiKA 



+^p(A^ - IjmgA'^K - 8vd{N - l)qyoX^K (7.30) 



and 

dt 



-2vd{N - l)/fy3't + ^Vd{N - l)/^y2«A 

+ lQvd{N - l)4yiKU3 - ^^{N - l)m'i {ul + u^X^ k 

+l^(Ar - l)miu^\^K + Avd[N - l)liu^yon 
-2Qvd{N - l)liyiK\^ - 2^Vd{N - l)liu^yQK\ 

-^^(A^ - l)m^A^K + 2^Vd{N - l)liyQK\\ (7.31) 
We recognize that (t/q = k'^~'^Z^'^Y) 

Zn = K-Vn (7.32) 

holds in leading order and that the system above exhibits fixpoints of the orders 

(y„),~A"+^ (7.33) 
These results agree with our assumptions made above. 

The only point missing for a proof that eq. (|7.7] ) is indeed the leading contribution for large A 
concerns the neglection of the Q^-dependence of Z}.{p, Q"^) and rj. If this momentum dependence 
contributes only in nonleading order in A one obtains 7 = — 1 for d = 4. It is also conceivable, 
however, that the momentum dependence modifies the value of 7 and this is the reason why 
we have kept 7 as a free parameter at the present stage. In principle, these questions can be 
answered by establishing an appropriate fixpoint behaviour but we have not performed this 
task in the present paper. It may even be hoped that rigorous bounds on the maximal size of 



the various quantities appearing on the r.h.s. of (7J) could be proven with the use of exact 
evolution equations, in a spirit similar to our discussion for U3, z\ and z^- Even without going so 
far we believe that our estimate ( |7.7| ) for (i\ is reliable for very large values of A, with a certain 
uncertainity in the region of intermediate A. (There is disagreement about the behaviour of /3a for 
large A in the literature |18|. These references find the asymptotic behaviour /3a ~ A, /3a ~ A^ and 



/3a ~ A2 .) Excluding ji to be almost zero for intermediate A this behaviour is sufficient to derive 



the "triviality bound" ( 7.16 ) for the mass of the Higgs scalar. We emphasize that our treatment 
establishes the existence of a triviality bound without the use of a lattice regularization. For 
example, we could use a smooth ultraviolett cutoff A within our formulation in continuous space. 
This determines the "initial value" Fa for a given bare theory. Triviality bounds appear for a 
wide class of bare actions for wich Fa obeys ( [7.26D and ( [7. 27)) and similar restrictions for the 
momentum dependence of the vertices in Fa- 

8 Conclusion 

We have computed the /3-function in order ~ A^ and the anomalous dimension in order ~ A^ 
without ever doing a two loop calculation. We used instead an exact evolution equation which 
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has a one loop form, with classical propagator replaced by the full propagator and solved it 
using a nonperturbative method which combines exact evolution equations for "independent" 
couplings with renormalization group improved one loop expressions for "secondary" ones. It 
is impressive to see how different properties of the propagator in a background field combine to 
reproduce the perturbative two loop result in next to leading order in the small quartic coupling. 
We emphasize that in the course of our calculation we have in addition gained a lot of infor- 
mation on the momentum dependence of various IPI Green functions. In contrast to the more 
formal perturbative loop expansion the different contributions to the anomalous dimension and 
the /3-function for the quartic coupling are here directly connected with physical properties, i.e. 
the behaviour of Green functions. 

Even though the reproduction of perturbative two loop results by a one loop calculation is in- 
teresting by its own, this was not the main motivation of the present work. The two loop results 
are known since a long time and the present method does not offer calculational simplification 
in this respect. The main interest lies in the fact that our method is by no means restricted to 
an expansion in a small coupling A. The exact evolution equation is valid for arbitrary values 
of the quartic coupling. We have demonstrated this by computing the /^-function also for very 
large values of A. Our results confirm the lattice results for the triviality bound for the Higgs 
scalar mass in the context of a Euclidean field theory formulated in continuous space. 
A computational method giving reliable results also for large A is particulary important in di- 
mensions less than four, where it is well known that the critical behaviour is governed by a large 
value of the infrared fixpoint for A. Our method can be directly applied to arbitrary dimensions 
and arbitrary values of the coupling. The main modification as compared to the present calcu- 
lation is the size of the dimensionless mass term 2\k in the propagator. It cannot be treated as 
a small quantity anymore and the threshold functions have to be evaluated numerically. In view 
of the nonperturbative applications it is crucial to have a consistent nonperturbative scheme for 
an approximative solution of the exact evolution equation. This was presented in the present 
paper and amounts to a renormalization group improved one loop calculation of the " secondary" 
couplings. The successful precision test of our systematic truncation scheme offers the exciting 
prospect of precision calculations for the critical behaviour in two and three dimensions or the 
high temperature phase transition in four dimensions. 
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A One loop calculation of the wave function renormalization 

In this section we perform a one loop calculation of the wave function renormalization. From 
eq. ( |3.1| ) we get after functional derivation and inserting a configuration with constant scalar 
field 991 = \f2p, = for a 7^ 1 



{vf{Q-Q-p) 



ah 



U'k{p)Sab + 2pU'^{p)Sal6bl 



(A.l) 



From this we see 



and 



Q^Zk{p,Q^) 
Q^Zk{p,Q^) 

Zk{p,Q^ 
Zk{p,Q'^ 



{r^^\Q,-Q;p))^^-{r^^\o,0;p) 

(r(^)(Q,-Q;p))^^-(ri^)(0,0;p) 



22 

11 



(A.2) 



0) 
0) 



lim — — ?r 



(rl2)(Q,-Q;p))^^, 



(A.3) 



We consider the difference 
,2. ^ 1 



Zk{p,Q')-Zk{0,0) 



(rf (Q,-Q;p))^^-(rf^(o,0;p) 
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{r'i\Q,-Q;0) 
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(A.4) 



as one of the "secondary" couplings F^"'''^^^ in the sense of section 2 and similar for Z^- (Note 
Zk{0, 0) = Zk{0, 0).) In order to evaluate the one loop formula (|2.6|) we approximate F^^'* on the 
r.h.s. by the second functional derivative of the expression 



This yields 

Zk{p,Q^) 
Zk{p,Q') 



Zk{0,0) 



FfcM = J d^'xI^Ukip) + ^ZkO^if'^d'^^a 

Uf{p)p 



(A.5) 



{27r)-'J d\ {P{q) + U'kip) + 2U'^{p)py' 



{P{q + Q)+ U'kip))-' - (PiQ) + U'k{p)r' 



Zfc(0,0)-^(2^) 



{W'^{p) + 2U^'{p)pY 
P{q) + U',{p) + 2U'^{p)p 



{P{q + Q)+ U'kip) + 2U'^{p)p)-' - (Piq) + U'kip) + 2U'^{p)p)-' 



+ 



{N-l)Uf{p) 
P{q)+Ul{p) 



[P{q + Q)+Ul{p)r' -{P{q)+U',{p)) 



-1 



(A.6) 



and 



Zk{p,Q) 



dxxip^P + U',) \p + Ul + 2U'^p)- 



Zk{0,0) + 8^ufp 

Zfc(0,0)+4^(3f/^ + 2O)V /„ dxx^2P^P + U'^ + 2U'^p)- 



+4(iV 



dxx2P^ {P + U'k] 



(A.7) 
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where x = q ,P = ^ and P is given by (|3.25D . The one loop expressions eqs. (|A.6D and eqs. 
( A. 71 ) are ultraviolet finite in any dimension less than six. The renormalization group improved 
one loop approximation should therefore only be used for d < 6. 



B Evaluation of some integrals 



In this section we evaluate the integrals /(Q^) and J{Q^) appearing in section 5 in four dimen- 
sions. Throughout this section we will use the inverse propagator ( p.25|) for computations, i.e. 
we will set Z = 1. We start evaluating the I{Q^) which is defined in eq. ( p. 81 ) and reads, with 
a = k'^a, b = in eq. ds] 



HQ') 



1 



{27rr^ d\p-Hq) p-^{q + Q)-p-\q) 



1 1 



dadb e - 1 (a + 5) 



\-2 







(B.l) 



We transform y = -^j^ , dy 



—da and perform the y-integration 



m') = VA 



dbb- 



l — e "P" 1+6 



fc2 



6 + i;4ln2. 



(B.2) 



In terms of z = one finds 



I{Q'') = I dz Iz-^ - ^z-^ - z-^e~^ 



(B.3) 



Writing the integrand as a Taylor series one finally obtains 



-1 



\n+l 



^^(n + 2)!(n + l) V2 



1\' 



n+l 



fc2 



n+l 



which is eq. ( |5.1[1| ). Analogously we proceed with the J((5^)-integral. It reads 



(B.4) 



J(Q2) ^ -{27:)-^ J d\P~\q)[p-\q + Q)-P~Hq) 



V/^k / dc 



ln(2 + c) - 2 ln(l + c) + Inc + 



l+£ 

2+c 



dzz 



1 — e fc^ 



l + c 



. (B.5) 



We consider the double integral K and integrate by parts 



l + c 
1 2+c 



K = 



dc I dzz ^ ( 1 



In 2 + dc 



-1 



Cf c \ a 2 + C l + c Q 



l + c 

2+c 



g fc2 l + c _|_ 



gfc2 2+c _(_ 

l + c F 



dz „q! 
-e fc^ 



l+c 
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The brackets [• • •] contain four terms. In the first one we substitute z = and expand in a 
Taylor series 



Similarly we proceed with the next two terms 







2 



2 

In the last term (the double integral) we expand the integrand in a Taylor series and perform 

l + c 
2+c 

the / dz-integration for every term. Again we substitute z = and z = 







l+c 



Summing up these contributions, one finally obtains 

j(g2) = g ^ — Lil! — ^ (9lY^' I dz-^^ - , ^""'l, 

^Jn + 2)!(n + l) l^fcsy (1 - Ji {I - z)^ 

which is eq. (|5.1l|) . 



(B.9) 



' dz-^ } (B.IO) 



C Contributions to (3x 

Finally we add up all the contributions wich enter in the /3;^-function ( [3.23 ) in order A'^. Inserting 
eqs. (|1|), ( |6T7D and (|6]2|) in { ^1^ we write 



'I, 



2v4{N + 8)X^ - 4vlX^{ {WN + 44) 

(-ir 



(-ir 



2;^g(n + 2)!(n + l) V2 



-1 



+2E 



6+2n 



^Q(n + 2)!(n + l) 
-(iV + 2) (Sm^/f-^ 
From eqs. ( ^l0|) and (|3^ we find 



dz 



n+l 



{l-zf 



2 



n+l 



(1-z)^ 



(C.l) 



li=n dyy-^e-y (l - e"^) 



n-l 



(C.2) 
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where y = x/k'^ and from eq. (|1.3D . These integrals may be integrated by parts and lead to 
F-functions. From this we get 



1"^ 

n 


- Z"' - 1 




= ''4 




= (n + 2)!, 


jQ+2n 


= 2(n + l) 



'-^^-T^. (C.3) 

2"+! ^ ' 

From eqs. and (|6l5|) we obtain for n > ^ + 1 with dt {fiP'^) = -2xd^ {j^P^) and 
integration by parts 



d 



5„,|+i + (n - 1 - -)ei 



dp 
dy 



with 



piy) = TJi 



y 



Directly we find 

and after further integration by parts 



/c2 i-e~y' 

mg = 1 

. 1 



(C.4) 

(C.5) 
(C.6) 
(C.7) 



(Note that the ej^ are closely related to the m^.) Inserting these results in eq. ( |C.lJ ) it remains 
performing two sums, namely 



(-1)" 



2,(n+l)(n + 2)- 



(-1) 



^ n + 1 ^ 



n + 2 



with X = ^ , X = I and 



n=0 



dz 



(1 - zf 



n=0 ' n=0 

ln(l + x) + [lii(l + x) — x] 



(C.8) 



ln(l + z) 1 



1 - z 



2ln(l + z) + -ln(l-z^ 



(C.9) 



Putting all this together leads to the desired result. 
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Figure Caption 

Fig. 1: Graphical representation of the exact evolution equation 
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